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ABSTRACT 

We derive the asymptotic mass profile near the collapse center of an initial spherical 
density perturbation, S oc M"*^, of collision-less particles with non-radial motions. 
We show that angular momenta introduced at the initial time do not affect the mass 
profile. Alternatively, we consider a scheme in which a particle moves on a radial 
orbit until it reaches its turnaround radius, r*. At turnaround the particle acquires 
an angular momentum L = C^/GNL/rl per unit mass, where is the mass interior 
to r*. In this scheme, the mass profile is M (x r^/('^+3'^) for all e > 0, in the region 
r/rt <^ C, where vt is the current turnaround radius. If £ <C f then the profile in the 
region C<^r/rt<^isMcxr for e < 2/3. The derivation relies on a general property 
of non-radial orbits which is that ratio of the pericenter to apocenter is constant in a 
force field fc(i)r" with k{t) varying adiabatically. 
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1 INTRODUCTION 

In the hierarchical scenario for structure formation, the ini- 
tial density field is Gaussian with a fluctuation amplitude 
that decreases with scale (Peebles 1980). Non- linear gravita- 
tional evolution then causes matter to aggregate into bound 
virialized objects (halos) which are believed to harbor galax- 
ies, and galaxy groups and clusters. The mass distribution in 
halos can be inferred from a variety of observations. These 
include observations of rotation curves of spiral galaxies 
(e.g., Persic et. al. 1996), motions of satellite galaxies (Zarit- 
sky & White 1994) , lensing distortions of background galaxy 
images by the potential wells of massive halos (e.g., Bartel- 
mann & Schneider 1999) , and velocity dispersion and X-ray 
maps of galaxy clusters (e.g., Bahcall & Fan 1998, Sarazin 
1986). Nevertheless, no single observational method probes 
the mass profile over the entire extent of the halo. Therefore 
a conclusive analysis of the observations must rely on an 
assumed model for the halo mass profile. This is important 
for example in estimating the total masses of galaxy clus- 
ters. Further, halo profiles in the inner regions depend on 
the type of dark matter (e.g., Moore 1994). A comparison 
between observed and model mass profiles in the context of 
a given cosmological model can provide important informa- 
tion on the nature of dark matter. 

We lack a complete theory of nonlinear gravitating sys- 
tems. So a detailed study of nonlinear collapse into bound 
objects must rely on N-body simulations. The simulations 
seem to be converging on the shape of the mass profile in the 
halo inner regions (Navarro, Frenk & White 1997, Moore et. 
al. 1998, Jing 1999, Klypin et. al. 2000). As of yet no satisfac- 



tory analytic explanation of the simulations results has been 
suggested (but see Weinberg 2000, Syer & White 1998). In 
a generic collapse the infalling matter is in irregular clumps 
undergoing dynamical friction and tidal stripping as they 
sink towards the center (e.g., Avila-Rees et. al. 1998, Nusser 
& Sheth 1999). Moreover, a collapsed object today might 
have gone through a merger with another object of compa- 
rable mass, a process that might have an effect on its current 
mass profile. Because of these factors, the focus of analytic 
studies has been the collapse of perturbations with special 
configurations. An important step forward in these studies 
came with the realization that the potential in the inner 
collapse regions varies very little over the dynamical time 
scale (Gunn 1977), admitting an adiabatic invariant for the 
motion of a particle. Using adiabatic invariance, Fillmore 
& Goldreich (1984, hereafter FG84, see also Bertschinger 
1985) have shown that in spherical symmetry an initial den- 
sity perturbation AM/M ~ develops an asymptotic 
mass profile M oc r near the origin for < e < 2/3 and 
M cx r-^/(^+^") for e > 2/3. FG84 restricted their analy- 
sis to particles collapsing on purely radial orbits. Spherical 
collapse with non-radial particle motions has also been con- 
sidered in the context of adiabatic invariance. Gurevich & 
Zybin (1988a, b) developed a formalism based on the the- 
ory of adiabatic capture (e.g., Lifshitz & Pitaevskii 1981), 
to estimate the mass profile in a spherical collapse including 
non-radial motions. Although their formalism can be useful 
to study the collapse of various initial conditions, Gurevich 
& Zybin applied it only to the collapse near the center of 
a smoothed initial density maxima. Ryden & Gunn (1987) 
included non-radial motions in their application of adiabatic 
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invariance in a numerical scheme to study the evolution of 
a spherical peak in a Gaussian density field (Bardeen et. al. 
1986, Hoffman & Shaham 1985). 

White & Zartizky (1991) described the M cx r profile 
for e < 2/3 in the purely radial spherical collapse in terms of 
the crowding of orbits as particles pass through the center. 
They conjectured that if particles had non-radial orbits then 
the crowding effect is avoided and the scaling M oc r^^^^+^'^'i 
e > 2/3 would also be valid for < e < 2/3. If ere we examine 
this conjecture in detail by generalizing the analysis of FG84. 
We consider two schemes for assigning angular momenta to 
the collapsing particles. 

The outline of this paper is as follows. In section 2 we 
write the equations of motion and the relevant initial con- 
ditions and present two schemes for assigning angular mo- 
menta to particles. In section 3 we discuss self-similarity and 
adiabatic invariance in the context of non-radial motions. In 
section 4 we derive the asymptotic profiles for various cases. 
We conclude with a summary and discussion in section 5. 



2 THE EQUATIONS 

We consider the evolution of an isolated spherical positive 
density perturbation in a flat universe made of collision-less 
matter only (density parameter SI = 1). We describe the per- 
turbation by a large number of equal mass particles having 
angular momenta in random directions such that the mean 
angular momentum at any point in space is zero. Spherical 
symmetry is then preserved and the angular momentum of 
each particle is conserved as the system evolves. It is often 
convenient here to think of a particle spherical shell 

obeying the equations of motion of a point particle moving in 
the gravity fleld of the perturbation with conserved angular 
momentum. We write now the equations of motion govern- 
ing the evolution of a shell with a given angular momentum 
per unit mass, L. Denote by r be the distance from the sym- 
metry center, and by M(< r,t) the mass contained inside r 
at time t. The trajectory r{t) of the shell as a function of 
time is then given by the solution to 

GM{<r,t) 



dV 
dt2 



-I- 



(1) 



where the mass M{< r,t) is determined from the distribu- 
tion of shells at time t. Deferring the discussion of the way 
shells (particles) are assigned angular momenta, these equa- 
tions determine the evolution of the system given a set of 
initial conditions specified at an early time ti — > 0. We as- 
sume that the initial radial velocity of a shell at Vi is equal 
to the Hubble expansion velocity at ti, that is, 

^ - -r (2) 

We express the initial mass distribution in terms of the rel- 
ative mass excess Si = AM/M interior to an initial radius 
ri = [SM/A-Kptiti)]^^^ where pt{ti) = (evrGt)"^ is the mean 
background density ti. We take the scale free form 



AM _ /My 
IT ^ VMo/ 



(3) 



where e > 0. A perturbation with e > 1 can be realized 
by placing a point mass at the center of a void with local 
density contrast cx — M"" (cf. Chuzhoy & Nusser 2000 ). In 



cosmology the initial perturbations are small so we restrict 
the analysis to shell far enough from the center so that Si <C 
1. 

The general solution to the equations can be obtained 
by numerical integration where the mass distribution is con- 
tinuously updated from the new positions of the shells at 
each time step. However the equations allow an analytic so- 
lution for the motion of a shell that has not crossed any 
other so that M[< r{t)] — const. We will see later that for 
some choices of the angular momenta, this solution is appli- 
cable in the outer regions of the collapse. The solution can 
be written in the parametric form (e.g.. Landau & Lifshitz 
1976), 



r = ro(l — ecosT^); 
where 



1 + 



2El/_ 



ro 



t = to{ri — esiuT^) 



to 



GM 1 



GM' 



(4) 



(5) 



According to this solution, the shell expands until it reaches 
a maximum expansion radius (turnaround radius) r, — 
ro(l -I- e), at time t, = nto. The solution ceases to be valid 
shortly after maximum expansion when the shell crosses 
other shells returning from the inner regions after passing 
through their maximum expansion sometime ago. The ex- 
pressions for r* and for purely radial motion can be re- 
covered by taking the limit L ^ in (^. In this limit we 
find (cf. Peebles 1980, FG84) 



ro{l+e) = 



GM 



GM(1 
riS~^ = 



e) 

3Mo 



1 1/3 



GM 



3/2 _ 



^U5 



- ^ 
4 

37r 



'iTVpb{ti 
1 



(-) 



1/3-l-e 



67vGpb(ti) 
M \ 3^/2 



1/2 



-3/2 



37r fM_\ 



(6) 



(7) 



(8) 



(9) 



The solution (^ also describes the motion of a parti- 
cle with non- vanishing angular momentum in an attractive 
force field oc 1/r^. We remark here that by taking the an- 
gular momentum to zero we do not recover the motion of a 
particle with zero angular momentum. In the limit L ^ 
the particle has zero radial velocity near the center oscil- 
lates between r = and a maximal radius, r,. A particle 
with L = on the other hand has an infinite radial velocity 
at r = and oscillates between — r, and r, . 

2.1 Two schemes for assigning angular momentum 

We will consider two schemes for assigning angular momenta 
to particles. In the first scheme (scheme A) angular mo- 
menta are assigned at the initial time in such a way that no 
additional scale is introduced in the collapse. The angular 
momentum of each particles is conserved in the subsequent 
evolution of the perturbation. In the second scheme (scheme 
B) a particle acquires its angular momentum when it is at 
maximum expansion, prior to that it is assumed to have a 
purely radial motion . The angular momentum of a particle 
in scheme B is oc per unit mass (White & 

Zartizky 1991), so no additional physical scale is introduced. 
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In scheme A the energy per unit mass of a sheU can 
be written in terms of the initial radius and mass excess 
SiMi inside the shell at the initial time, ti, 



2r? 



G 



S^Mi 



(10) 



so to make all energy components scale similarly we choose 

1 1/3 



= 2aGS^Miri = 2aG 



Mf'sf'-^ , (11) 



where 1 > a > insuring that the energy E = [a ~ 
l)G5Mi/ri is negative and the shell is bound. This choice 
for L does not introduce any scale in the initial conditions. 

The eccentricity, e, corresponding to the motion before 
shell crossing is, according to given by 



e = 1 + 4a(a - 1)5' 



(12) 



Since in cosmological perturbations the initial density con- 
trast is small, the last relation implies that the eccentric- 
ity is very close to unity. Therefore, according to (^), the 
turnaround time of the shell is = tiS~^^^ ~ rf'^ and the 
turnaround radius is r, = rid~^ ~ '''l^^" ■ This scaling means 
that inner shells collapse earlier than outer shells and the so- 
lution ^ is always valid outside the the radius of the shell 
at maximum expansion at the current time, t (hereafter, the 
current turnaround radius). 

In scheme B, a shell moves with zero angular momen- 
tum until it reaches its turnaround radius, where it is as- 
signed an angular momentum per units mass according to 



'lnpb{ti 



1/3 



^4/3^4/3+. 



(13) 



where £ is a constant value for all shells. To insure that a 
shell remains bound and subsequently collapses towards the 
center we demand £, < 1. Note that the angular momentum 
assigned to a particle at the initial time in scheme A is oc 

Shortly after reaching its turnaround radius, a shell 
starts its oscillations in the collective gravitational potential 
well generated by shells that have passed their turnaround 
radii at earlier times. Angular momentum introduces an ef- 
fective repellent force, L^/r^, preventing shells from reach- 
ing the center. We will see that in scheme A angular mo- 
mentum affects the density profile only inside a radius equal 
to the initial radius of the shell at the current turnaround 
radius. For cosmological initial conditions this radius tends 
to zero since ti 0. So angular momentum does not play a 
role in fixing the density profile of the evolved perturbation. 
In this case the results of FG84 for collapse of particles on 
purely radial orbits remain valid. In scheme B, angular mo- 
mentum can prevent particles from penetrating a significant 
fraction, depending on C, of the current turnaround radius. 



which at any time t can be taken as the current turnaround 
radius, rt{t). The mass scale corresponding to the current 
turnaround radius is Mt = M{< rt). So the mass distribu- 
tion M{r,t) must satisfy the self-similarity condition 



M{r,t) = Mt M{r/rt) 



(14) 



where TVf is a function of r/rt only. Substituting tt, — t in 
(0) we find 



rt 



3Mo 



1 1/3 



^ N 2/3+2/(9e) 



Mt = M, 



t \ 2/(3£) 



(a 



(15) 



(16) 



Assume that in some region the mass distribution can 
be approximated as 

M = kity, k{t) = kor" , (17) 

where ko,s, and 7 are constants^ We will determine 7 in 
the next section by generalizing the analysis of FG84 for 
collapse with no angular momentum (see also Zaroubi & 
Hoffman 1992). However, the following useful constraint 



s + 



3e 



7 1 + 



3e 



: 



(18) 



can readily be found by substituting the asymptotic form 
( p^ for the mass in the self-similarity condition (^^ and 
using (^. 

Throughout the rest of the paper we will assume that 
the potential well near the center varies adiabatically (Gunn 
1977, FG84). This means that a shell near the center makes 
many oscillations before the potential changes significantly 
(Gunn 1977, FG84). In a slowly varying potential the action 
variables associated with the motion of a particle is invari- 
ant. The angular action variable is the angular momentum 
which, thanks to spherical symmetry, is conserved indepen- 
dent of whether or not the potential changes adiabatically. 
The radial action variable associated with the motion of a 
particle in the potential well of (^^ is 



27r 



dr 



dr 



2 \ E- 



7-1 



1/2 



(19) 



(20) 



where and rt are, respectively, the apocenter and pericen- 
ter of the particle's orbit. We prove now that the invariance 
of Jr implies that the ratio ^ = rt/ra is constant with time. 
This important property will simplify the determination of 
7 in the next section. We prove it in the following way. The 
energy E and k can be expressed in terms of ra and rt as 



2k 



7-1 ' 



(21) 



3 SELF-SIMILARITY AND ADIABATIC 
INVARIANCE 

The initial conditions are scale free, the two schemes for as- 
signing angular momenta do not introduce any additional 
scale in the problem. So the only characteristic scale in the 
evolution of the perturbation is the scale of non-linearity 



2E = -L^ 



Therefore 



.7-1 



' a 



(22) 



* For collapse described by a finite number of particles, the form 
JitI ) is meant to describe the mass averaged over a few particle 
oscillations within r. 
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Jr = 2tvL 



du 



(r^ - 1) 



,7-1 



+ (1- 



1/2 



(23) 



The angular momentum L is conserved, so the invariance of 
Jr means that ^ is constant. Since the only special time in 
the particle history is its turnaround time , the invariance 
of ^ implies that 



ra,b 



(24) 



near the center. The index q can be expressed in terms of 7 
and s by using ( p^ ) in ( ^ ) . This yields 



g(7 + l). 



(25) 



4 THE ASYMPTOTIC BEHAVIOR 

Assuming the asymptotic form (^) we now turn to es- 
timating the exponents 7, s, and q. Let P{r,ri,t) = 
° dr( dr/ dt)~^/ JJ^" dr( dr/ df)~^ the fraction of time a 
particle with initial radius n spends inside the radius r at 
the present time. The function P is well defined if the mass 
within r changes very little during one oscillation period of 
a particle inside r. The mass M{< r,t) can be written in 
terms of P as 

M{< r,t)= / dM,P{r,r„t) , (26) 



where Mi = 4npt{ti)rf /3. By substituting (0) and (0) in 
this last relation we obtain 



dM. 
Mt 



P{r, n) 



(27) 



Following FG84 we write P in terms of u = r/va- By def- 
inition, P{u) — for u < ^ and P{u) = 1 for m > 1. For 
^ < u < 1, we write P{u) — /(it)//(l), where 



/(n) 



dv 



(r 



.,7-1 



1 



1 -1/2 



1 



Using the relations and (^ 
terms of u. So 



^7-1 

I we express the mass, Mi, in 



(28) 



where 



duu 



-(i+p) 



r/rt 



(29) 



(30) 



^ 2 + 3(2-3g)e ' 

Particles with rt > r, i.e., u — r/ra < ^, have P{u) = and 
do not contribute to the integral on the right hand side. The 
relation (|2^) allows us to obtain the asymptotic behavior in 
scheme A and B without a detailed calculation of P(it, ,^). 

4.1 Scheme A 

The first step in the derivation of the asymptotic behavior 
is to relate ^ to the initial density contrast. This can easily 
be done by substituting the expression (|l^) for L in (pl|), 
and using (F3) and (p|). The result is 



7+1 



e-i 



-a(a-l)(7-l)5? 



(31) 



Since Si <ti 1 this relation implies that ^ ^ 1 as well. Re- 
taining only the lowest order in 5, the relation gives ^ oc Si 
for 7 > 1, and ^ oc s'^^^~'^^'> for 7 < 1. Since < r», we 
have rt — < £,r, — ^ri5~^ < ri, i.e., the pericenter 
of a shell is smaller than its initial radius. Except near rt, 
where a shell spends only a tiny fraction of its orbital time^ 
the motion is almost radial. The initial radius, ru, of the 
shell currently at turnaround sets an upper limit on rt- For 
r ^ ru all shells move on almost purely radial orbits and 
the analysis of FG84 for collapse on purely radial motions 
apply. Therefore, according to FG84, we have 



7 = p = 



and, 



7=1 



1 + 3e 



s = g = 0, 



for 



- 3 ' 



3e 



2g, 



for 



2 

^<3 



(32) 



(33) 



However ru tends to zero as ti 0. So for cosmological 
initial conditions in which ti —> the the collapse is identical 
to that of particles on purely radial motions. 

4.2 Scheme B 

The angular momentum in this scheme is given by (|l|). So 
the relation (|2l[) gives 



7+1 



= -£^(7-1) 



(34) 



implying that ^ is the same for all particles. The radius 
rt(, is an upper limit on the pericenters of all shells. So 
only particles which have passed their turnaround radii early 
enough can contribute to the density at r < rt^. Particles 
at r ^ rj^^ have their apocenters inside rt^ as well, and 
so the radius rt^^ marks the boundary of the inner region 
where we expect the same power law behavior for the mass 
distribution. In this region only particles with rt = ^r^ < r, 
i.e., u> ^ contribute to the integral in (pa), so 



r\-<-P 1 



(35) 



The integral is independent of r, therefore we must have 
7 = p for all e > 0. Using the expression (js^), and the 
relations (|l|) and (^, yields 

3 

"/ — p = , s = q = , for all e > . (36) 

1 + 3e 

In the collapse with L — (FG84) and in scheme A this 
result is correct only for e > 2/3. 

If 5 <C 1 then in the region rt ^ r 3> rt^ particles 
move on almost radial orbits and the analysis of FG84 for 
collapse with no angular momentum is applicable. Applying 
the analysis of FG84 in that region gives 7 = 1 for e < 
2/3 and 7 = p, as in (p6|), for e > 2/3. Therefore, as r 
is increased, the mass profile index varies from 7 = p to 



t Consider the orbit of a particle in the plane (r, ip) where ip is 
the angular position. The time spent near ri, can be estimated as 
< 27rr'^/L, where is a few times r^. By 

Jill ) we find that r^/L ~ ^"Tc where n > and Tc is the period 
of a circular orbit at the apocenter ra- So the time fraction spent 
near is negligible for ^ <C 1. 



© 0000 RAS, MNRAS 000, 000-000 



Self-similar spherical collapse with non-radial motions 5 




Figure 1. Mass profiles from spherical N-body simulations in 
which particles acquire angular momenta according to scheme 
B. The initial density perturbation has e = 0.3. The solid, long 
dashed, short dashed, and dotted curves are, respectively, mass 
profiles for £ = 0.9, 0.3, 0.1, 0.001. 



7 = 1 for e < 2/3, but remains 7 = p for e > 2/3. In 
figure 1 we show the mass profiles obtained from spherical 
N-body simulations of 12000 shells for e = 0.3 with angular 
momentum introduced according to scheme B. The curves 
correspond, respectively, to simulations with C = 0.9 (solid), 
0.3 (long dashed), 0.1 (short dashed), and 0.001 (dotted). To 
save CPU, a shell participates in the simulation only after it 
reaches its turnaround radius. A fourth order Runge-Kutta 
time integration is used to integrate the equations of motion. 
The mass profiles in the figure behave as predicted from 
our analysis. A particle with £ — 0.9 moves in an almost 
circular, so that it effectively circles the halo at a radius 
slightly smaller than its turnaround radius. Therefore the 
corresponding mass profile (solid curve) would have an index 
close to 7 = 3/(1 + 3e) = 1.57 almost for all r < vt- The 
mass profile in the inner regions for £ = 0.3 and 0.1 (long 
and short dashed, respectively) have 7 ~ 1.57, while in the 
outer regions it becomes « 1. For £ = 0.001, particles move 
on almost radial orbits over the whole distance range shown 
in the plot. Hence 7 ~ 1. 



5 DISCUSSION 

We have studied the effect of non-radial motions on the mass 
profile of spherical self-similar collapse. Assigning angular 
momentum at the initial time does not affect the evolved 
mass profile. This is not surprising as the initial kinetic en- 
ergy is dominated by the radial Hubble expansion velocity 
oc r/ti. Any additional finite velocity is negligible compared 
to the Hubble expansion velocity as ti —> 0. This has also 
been shown explicitly for collapse with purely radial motions 
where a finite radial velocity component is superimposed on 
the Hubble expansion at the initial time (e.g., Peebles 1980, 
Padmanabhan 1993). We have also considered the conse- 



quences of a different scheme (scheme B) for assigning angu- 
lar momenta to particles. In scheme B, a particle acquires its 
angular momentum at maximum expansion. The acquired 
angular momentum is chosen in such a way that the ratio of 
the pericenter to apocenter is the same for all particles and 
no additional physical scale is introduced. In this case we 
have confirmed the conjecture by White & Zaritsky (1991) 



that the scaling M oc r' 



,3/(l-|-3e) 



is maintained for all e as 



r ^ 0. But we have also shown that if the added angular 
momentum is small, then the behavior M oc r is restored 
far enough from the center for a perturbation with e < 2/3. 

What is the justification for scheme B in a generic col- 
lapse configuration? Consider a halo forming in a gravita- 
tional collapse of a high density region in an initial scale 
free random Gaussian. The infalling matter is distributed 
in clumps (satellites). If the halo is approximately spherical 
then the collapse process cannot introduce any additional 
scale and angular momentum must be proportional to the 
product of the mass and the radius at the epoch of maximum 
expansion. But scheme B also assumes that the angular mo- 
mentum of a particle is negligible before it reaches maximum 
expansion. This is a reasonable assumption if angular mo- 
mentum is mainly generated by tidal interaction between the 
infalling satellites themselves. Because the halo is assumed 
spherical, the halo-satellite tidal interaction is not expected 
to generate any angular momentum with respect to the halo 
center. The assumptions underlying scheme B are consis- 
tent with results of N-body simulations (e.g., Barnes & Efs- 
tathiou 1987). Adiabatic invariance was also applied to non 
scale free initial density perturbations (Hoffman & Shaham 
1985, Ryden & Gunn 1987, Lokas 2000, Lokas & Hoffman 
2000). N-body results imply that in these situations as well, 
it is reasonable to introduce angular momentum according 
to scheme B. 
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